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The Randall- Sundrum model is studied in 6 dimension with AdS4 or dS4 metric in the physical 
4 dimensional space. Two solutions are found, one with induced 5-dimensional gravity terms added 
to the induced cosmological constant terms. We study the graviton modes in both solutions by 
transforming the mass eigenvalue equation to a Schrodinger equation with a volcano potential. 
The spectrum of gravitational excitations depends on the input parameters of the theory, the six 
dimensional and the effective four-dimensional cosmological constants. The model gives a physically 
acceptable spectrum if the 4 dimensional cosmological constant is sufficiently small. 
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^ ■ I. INTRODUCTION 

U 

Cl^' Recent intensive study of string theory revived interest in gravitational models with more than 4 space-time di- 
, mensions. Models with extra (compactified or infinite) dimension, confining particles on three-branes 0,0, 0,0] have 
been in the forefront of research in particle physics in the past few years. These models have features that, unlike 
that of traditional string theories, may make the observation of extra dimensions possible in the near future. One of 
the most successful models proposed by Randall and Sundrum (RS) 3, 4] is based on an orbifold solution of the 5 
dimensional Einstein equation in AdS space, with a 3-brane fixed at the zero of the fifth coordinate, u (and possibly 
at another value of u, as well), around which the metric has a Z2 symmetry. This model may provide a solution to 
the hierarchy problem of particle physics. It confines gravitons to the neighborhood of a brane. 
^ The RS model with infinite extra dimensions [4] was extended to I? > 5 dimensions by Arkani-Hamed, Dimopoulos, 

(— I Dvali and Kaloper The model has intersecting {D — 2)-branes with gravity localized at the intersection. At long 
distances, along the intersection, Newton gravity is recovered. With an appropriately chosen curvature Newton's force 
is modified only at sub-millimiter distances. The cosmological constants were fine tuned to result in a Minkowski 
space at the intersection of the (D-2)-branes. 

The original RS model is required to satisfy a fine tuning condition if we demand that the 3-brane of our world 
■ has Minkowski metric. However, in our world the cosmological constant is probably non-zero and could have had a 
, large value in the past. Thus, it is of interest to investigate the generalization of the RS model when the physical 4- 
Cn| ' dimensional subspace is not Minkowski (M4), but AdS4 or dS4. Such a model has been investigated in 5 dimensions [1, 
.^-^ . 0, SQ. In horospheric coordinates, the space is infinite for AdSs, while it is finite for dS^. For distances much smaller 
' than the 4 dimensional curvature radius and much larger than the five dimensional curvature radius Newtonian gravity 
00 ' is reproduced. 

. The aim of the present paper is to extend the D — 6 case of @, 0, S to a warped 4 dimensional subspace. Just 
^ I like Ref. we require Z2 x Z2 orbifold symmetry in the two extra dimensions. In addition, for simplicity, we also 
. impose Z2 symmetry under the exchange of the two extra coordinates. In a previous work Chodos and Poppitz [lo| 
' discussed a three brane embedded in a six dimensional bulk using a different type of ansatz. Their solution also 
incorporated a nonvanishing cosmological constant on the physical three brane but the symmetry of their ansatz is 
AdS4 X E{2). 

Kaloper also found a D=6 solution with dS or AdS geometries on the physical 3-brane [ll| . His solution uses two 
4-branes intersecting at an angle different from the right angle. The model has Z2 x Z2 symmetry. It is asymmetric 
in the exchange of the two extra dimensions. Upon the analysis of small fluctuations no confined graviton mode was 
found in the case of AdS geometry on the physical 3-brane. 

In our analysis we have found two new solutions, both with two intersecting perpendicular 4-branes in the action. 
They are symmetric to the exchange of the two extra coordinates. These solutions will be discussed in the next section. 
In Section 3 we investigate small oscillations around the global solutions. In a subset of solutions that includes the 
ground state we are able to transform the mass eigenvalue equation to a Schrodinger form, with a now well-known 
volcano potential. We also find the spectrum of Kaluza-Klein modes. In Section 4 we conclude our paper. Two 
appendices contain some details of the calculations. 
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II. SOLUTION OF THE EINSTEIN EQUATION 

In 6 dimensions, when the 4 dimensional physical space is Minkowski, a solution can be found for the orbifold 
problem Q, just like in 5 dimensions, using conformal coordinates. One cannot however transform the coordinates 
to a horospheric coordinate system that has been used to solve the the 5 dimensional RS model without fine tuning 

aa. 

Unfortunately, in 6 (or higher) dimensions, when the 4-dimesional space is warped, the metric cannot be brought 
either to the horospheric, or to the conformal form. Each of these coordinate systems would only contain a single 
symmetric function of the two extra coordinates, while, as we will soon see, the metric components depend on at 
least one non-symmetric function of these variables. In fact, using appropriate coordinate transformations the most 
general Z2® Z2® Z2 symmetric ansatz for a 6-dimensional metric which is AdS (dS) in the 4-dimensional subspace 
(AdS4 or dS4) that, with a further gauge transformation, can be brought to the form 

ds^ — jABdx^dx^ = fl^ \g ^ydx^ dx'^ + g{u, v)du^ + g{v, u)dv'^'\ , (1) 

where 9^,^ is a 4 dimensional metric satisfying the Einstein equation with a four dimensional cosmological constant 
A, which may be positive (dS4 space), zero (Minkowski space), or negative (AdS4 space). Finally, is a symmetric 
function that can be brought to the form of the conformal factor for M4 space [H , 

^= [l + ai\u\ + \v\W 

With this choice, for AI4 space, the metric would reduce to the form of Ref. ^ and = 1. 

The only unknown function in the metric, g(u, v), has no symmetry property for the exchange of the two coordinates. 
As we will see below, the global solution of the Einstein equation for function g(u^ v) contains a single integration 
constant. This constant will be fixed later using the Israel junction conditions. 

A. Global solution of the Einstein equation 

The nonzero components of the Einstein tensor are fj^^ with ji ^ t,x, y, z, Quu, Qw, and Quv The four components, 
Qnn provide identical second order partial differential equations in variables u and v. It is easy to show that these 
equations follow, due to a Jacobi identity, from the rest of the equations. Thus, the bulk equations Quu ~ —guu^e 
and Quv = are sufficient to find a unique (up to an integration constant) solution. The equation Qy^ = —gy^Ag 
follows from the first equation using the replacement u v. These equations contain only first order derivatives of 
the unknown function g{u,v). 

To simplify notations we will use variables u and v rather than |u| and \v\. We may assume for the purpose of 
finding bulk (6-dimensional) solutions that u, w > 0. The role of Z2 x Z2 symmetry is important only when one solves 
the junction conditions that we will discuss after finding the global solution of the Einstein equation. 

The details of the solution are presented in Appendix A. We just recapitulate the results. The solution of the 
system of equations is given in the form 

9{u,v) ^ (3) 
ii(cr) + A{b, a) 



where we defined a — u + v and 6 — u — v. A((5, a) is an odd function of S. S{a) is given by 

Sia) = . 



2 1 

(4) 
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where C is a yet undetermined integration constant. 
The Taylor series of A{d, a) in 6 is 

00 

.4(<5,a) = ^<52'=+iafe(a), (5) 

where 



fc=0 



ao((T) = 5"(cr). 



(6) 
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and ak satisfies the recursion relation 

^ fc-i 

^ ' ^ ^ m— 

The recursion relation fT]) can be used to generate ak in arbitrary order. The first few are listed in Appendix A. 

Note now that g{u^ v) is unique, except for the choice of the scale a and the integration constant C. One of these 
constants can be fixed by setting the scale of variables u and v. The other will be fixed by the junction conditions to 
be discussed later. 

As usual, the scale is set by demanding that the metric is Minkowski at u = w = 0. As A{0, 0) = wc find the 
following value for the conformal scale from ^ 

Note that a solution may exist for both positive and negative cosmological constants, as long as ^ — ^ + C > 0. 



B. Israel junction conditions 



Due to the requirement of Z2 x Z2 symmetry the derivatives of the components of the metric tensor are not 
continuous at u = and v = 0. This discontinuity will generate junction terms in the Einstein tensor, localized on 
the u — and v — planes. These junction terms are: 













-5{u)guu 



8a{l + au) g^y{u,v)\y=o{^ + o.u)'' 
9(0, u) 
8a(l + au) 

9(0, u) 
8a(l + av) 

5(0, 



+ {u ^ v), 



(9) 



To satisfy the Einstein equations at u = and at w = two 4-brane contributions are required to cancel ([9]). Brane 
contributions come from the variation of w = and v = 4-brane terms such as a tension term (i.e. a five dimensional 
cosmological constant term) and possibly of dynamical terms. As we will see later the junction conditions will impose 
three u and v independent constraints on the constants of the theory. One of these can be used to fix the integration 
constant C. Then the brane contributions must in general contain two additional constants. One of these is the 
cosmological constant, A5. We choose the other constant as the 5-dimensional gravitational constant, G5. Thus, we 
propose to add the following terms to the action: 



1 



SttGs 



d'^xdux 



gR'^^^ + A5 / d'^xduxf^ + (u ^ u) 



(10) 



where i?'-^-' is the 5 dimensional Ricci scalar constructed from the induced metric on the t; = or m = brane. It 
is known that adding a four dimensional gravity term to the five dimensional action changes predictions for post- 
Newtonian effects, contradicting observations [13], OIlJl- The question, whether a similar problem arises about our 
solution with induced 5-dimensional gravity will be taken up in a future publication. 
The variation of pI7|) gives 



Qui + 9^iA5 = 5mm 



A5-A(l + au) + -— — + — (1 + auW— - 
.gu,0 2 r(w,0 



K^-2X{l + auf + 



g{u,Q) 



(11) 



There are similar contributions on the u = Q brane, obtained from (jlip by the exchange {u <-> v). One of the 
advantages of maintaining the u v symmetry is that one needs to deal with two junction conditions rather than 
four. 

It is easy to see that when we substitute our solution g{u, v) into ^ and pTjl multiplied by 5{v) the sum of these 
contributions does not vanish at all values of u. This alone does not mean that the 6-dimensional junction term cannot 
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be canceled by the brane terms of ((TO)) . The gauge of the junction contributions and of the brane contributions may 
be misahgned. In the next section we will how that this is indeed the case. 

Another way to satisfy the junction conditions would be to add further invariants to the action, such as Lovelock 
terms, but one would need an infinite number of such terms in the Lagrangian to have a chance to satisfy infinitely 
many constraints imposed by the junction conditions. 



C. Coordinate transformation 



As the six dimensional action is gauge invariant, we can apply a gauge transformation to its contribution to the 
junction condition. We require that the gauge transformations satisfy the following conditions: 

1. The (u <~* v) exchange symmetry is maintained. By this requirement we reduce the number of junction equations 
to be satisfied from four to two. 

2. We require that the transformed coordinates are such that u' and v' vanish when u and v vanish, respecitively. 
This constraint is required by the junction equations. It implies that the position of the branes is unchanged by 
the gauge transformation. 

3. We will maintain the scale of the coordinates at m t; = 0, by requiring that the components of the gauge 
transformed metric tensor and g^y also tend to one when u,v —^ 0. 

These requirements lead to the following form of coordinate transformations 

u u' ^ h{u, v) 

y V' ^ h{v,u), (12) 

where h{u,v) is a function, h{0,v) — 0, h^u{u,v)\u=v=o ~ 1- 

We will show in Appendix B. that the two junction conditions can be satisfied by an appropriate choice of h(u, v). 
While that result is important, the physical consequences of our solution can be explored without the detailed 
knowledge of the form of the function h(u, v) provided h{0, 0) = and h,u{u, v)\u=v=o = 1. 

When one writes down the junction conditions one must take care of the fact that the action is varied with respect 
to QABi rather than qab- In other words, certain linear combinations of ^ enter the junction conditions. Denoting 

these appropriate linear combinations by Ai?|^^ the junction conditions at ?; = are 

^^yR% + ^Rfs^iv) = 0. (13) 

One obtains two independent junction conditions, one for choosing A — B = ^ and one for A = B = u (see 
Appendix B). The junction contributions are functionals of two independent single variable functions, /i(u) — h{u, 0) 
and /2(u) — h^y{u,v)\y=Q, in an algebraic manner. The only constraints on these functions are /i(0) = and 
— /2(0) — 1. We must satisfy these three independent constraints. Once these constraints are satisfied one 
can always find two functions, fi{u), which solve the two junction conditions. The function /2(u) enters the junction 
conditions in a very simple manner: its third power is a multiplier of the 6-dimensional junction contributions. Taking 
the two junction conditions at u = provides two constraints on the constants of the model. Taking the ratio of the 
two junction conditions eliminates /2(u) from the resulting condition. Taking a derivative of the resulting equation 
with respect to u, at u — 0, provides a third constraint. The function h{u,v) enters each of the three constraints 
only through the combinations /i(0) = 0, /{(O) and /2(0), all fixed by our constraints on the gauge transformation. 
All other derivatives of the junction conditions taken at u contain unconstrained derivatives of h{u, v) that are 
determined by these very equations. 

It is clear now why, in general, we need the two terms in (jTU]) (cosmological constant and gravity). Recall that 
besides the coupling constants G5 and A5 we still have one undetermined integration constant, C , (or, alternatively 
the conformal scale, a, which is related to C through ([8]). The three equations obtained from the junction conditions 
are just sufficient to fix these three constants. In fact, one could replace the 5-dimensional scalar curvature term of the 
Lagrangian by another gauge invariant term, such as a Lovelock term, but the scalar curvature term is the simplest 
choice. 

The two junction conditions taken at m = provide the following equations 

= -^-7^{2A-6a'-A5), (15) 
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TABLE I: Summary of parameters characterizing the two solutions. 



Solution 


a 




c 


C?5 


As (As) 




S{a) 


Solution I ±^ 
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As 
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ZaGQ 
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-lOA- f As 


-6A 


—rvj^ + 1 + r 


Solution n zby 


1 llA 
100 


_ Aa 

20 


17A 
300 


CO 


a 
ttGq 








where we used the expression, obtained from our solution, 



ag,u{u, 0)|„=o = A - lOa^ - (16) 
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The set of equations (|T4|) and (|T5|) have two solutions. 
1. Subtracting the two equations from each other implies 



2 A Ae , . 



6 20 

or in other words, the integration constant C = 0. Then we also have 

The third constraint equation, to be discussed below will also determine the 5-dimensional gravitational constant, 
G5. 

2. Note that if we take the limit G5, A5 ^ 00 (no scalar curvature term on the brane) such that A5 — As/SvrGs is 
finite then ([14]) and (fT5| become identical, providing an expression for A5 

As = (19) 

The third constraint equations will fix the yet undetermined a (or, alternatively, fix C). 

Now for solution 1. the mere presence of the dynamical term on the w = and u = branes, combined with the 
fact that the metric tensors depend on |m| and respectively, requires the existence of a non- vanishing tension on 
the u = V — Q 3-branes, just like in the 5 dimensional Randall-Sundrum solution. We denote the tension by A4. We 
obtain 

A4 = -6A. (20) 

If A = (non- warped space) the contribution on the 3-brane vanishes and the solution reduces to that of 

Let us consider now the third and final constraint equation. As we indicated earlier it is obtained from the 
derivative of the ratio of the two junction conditions, eliminating /i,i,(u, u)|i,=o from the equation. To reduce the 
obtained equations to a constraint on our constants we need to calculate derivatives of g{u, v) aX u — v — Q. These 
can be determined cither using the explicit form of g{u, v) from ^ or directly from taking derivatives of (|14p and 
(|15p and setting u = v = 0. To simplify matters we write down the the third constraint separately for the first and 
second solution. For the first solution (G5 ^finite) we obtain, using (|17p 



1700A^ - 40A5(80 + 27A6) + 3(500A2 + 340AA6 + 57A^) = (21) 

For the second solution (G5 = 00), after substituting we obtain for the third constraint 

100a^-llA + 5A6==0 (22) 

These equations can be used to calculate the rest of the constants. All the constants of the theory are determined 
by the parameters of the 6 dimensional world, Gg and Ag and of the physical cosmological constant, A. The value of 
parameters found for the two solutions are given in Table 1, where we use the notations r = —\/Qa? and w = 1 + aa. 

We can see from Table 1 that in fact, taking into account the two possible signs for a, we have four independent 
solutions, rather than two. Note that at a special choice of the parameter r, r = —38/41 the combination 50A+ lOAg = 
and A5 = 0. Then the geometry on the intersecting branes is asymptotically Minkowski. 
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D. The range of variables u and v 



The range of variables u and v plays an important role in the investigation of small deformations of the above 
solutions. The linearized deformations determine the spectrum of gravitational excitations and the nature of the 
gravitational force among massive objects. 

The form ([3]) of the solution implies that zeros of S{a) ± A{5, a) (for both signs) constitute boundaries in the (w, v) 
space. Suppose the equation of such a boundary is on the curve 5 = 5[a). Then 5{a) satisfies 

[g{u, v)]~^ = S{a) + A{6{a),a) = 0. (23) 

This equation can easily be solved. Consider that whenever ((23)) is satisfied, unless 5 = 0, [g{v,uy\~^ ^ 0. Then 
the equation Q^y = 0, (|A2p of Appendix A implies that g,y{u,v) must also vanish. Then, differentiating (j23p with 
respect to a and comparing it with the substitution of (5 5{a) into the equation g^v{u, w) = we obtain 5' {a) = — 1. 
Integrating this equation we obtain 5[a) — —a + 2c, where c is a constant. Owing to the definition of 5 and cr, the 
boundary of the domain is determined by a constant value oiu, u — c. In other words, {g[c,v\)~^ — for all admissible 
values of v. This constant can be determined by taking the limit v c. At this point S = A{6,a) = and a = 2c. 
In other words, the equation 

S'(2c) = (24) 

provides the corresponding value of c. 

Notice now that when g^u = g[u,v] has a pole at u = c g^y = g[v,u\ also has a pole at v — c. Then, whenever 
has a positive solution the range of variables u and v is the square \v\ < c. 

Let us examine now the solution of (|24p for our four metric solutions. When a > the boundary r = separates 
domains with or without solutions. However, when a < the corresponding boundary value is r = rg, where tq — —1 
for Solution I and tq ~ —50/33 for Solution II. Owing to the identites r — — 6a^A and 1 — r/ro — a^20A6, these 
boundaries separate domains with different curvature: r > corresponds to AdS4 and r < to dS4 while r > tq 
corresponds to AdSe while r < tq corresponds to dSe. Thus, the boundaries are: 

1. Solution I, a > 

The boundary, c, satisfies 

-r(l + 2ca)2 + l + r = (25) 
where r = -A/Ga^ > 0. Then for r > (AdS4) the range of variables u and v is 

H,kl<^|\/^-ll- (26) 

For r < (dS4) the range of |u| and of |t)| is infinite. 

2. Solution II, a > 
([M)) has the form 

33 17 

- r(l + 2caf + 1 + — r + — r(l + 2ca)^ = (27) 
50 50 

No analytic solution of this fifth order equation exists but at small positive values of r the range of u and v is 
infinite. Increasing r we obtain a critical value, above which the range becomes finite, because (1271) has a positive 
solution. This critical value is obtained from the coincidence of the zero ([27l) with the zero of the derivative 
S{(7). This happens at 

- ■ - -2 33^ 17 . ^ 115.5, (28) 




■"min 50 80 '"mm 



where 



20X 1/3 

17 



1 + 2aCmin. (29) 



If r increases, starting from rmin than the boundary value of |u| and \v\ decreases from Cmin(''min) = 
(l/2a)[(20/17)i/3 - 1] = -^0.2095 to Cniin(oo) = ^0.0538 . For negative values of r ^ always has a posi- 
tive root, so the range of variables is finite. 
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3. Solution I, a < 

The largest possible physical range of w = 1 + aa is < w < 1. There is a solution of the equation S'(cr) — 
only if r < —1. This corresponds to dSe space. The solution of the S{a) = equation provides the following 
bounds for 

For r > —1 the range of u and f is bounded only by the condition I u| + |d I < l/|a| (i.e. finiteness of the conformal 
factor) . 

4. Solution II, a < 

At w = 1 S ^ I. Atw = OS=l + 33r/50. Therefore, if r > ro = -50/33 (AdSe) there is no solution of 
the S{a) = equation (it is easy to see that there is no minimum of S{a) lower than 0). However, there is a 
solution of the equation S{a) = for r < —50/33 because S'(w = 0) < 0. So w > Wmin > 0. These values of r 
correspond to dSg space (Ag > 0). Just like in the case a > there is no analytic solution for w < Wmin- 



III. SMALL OSCILLATIONS 



Suppose the metric solutions we previously found, denoted ^gab, are modified by a small perturbation, such as 

5a6(e) = °5afc + hab- (31) 

As now it is standard in the literature, we will impose the TT (traceless, transverse) gauge conditions on hab- In 
addition we will impose axial gauge conditions, hau = ^au = 0- Thus, the nonzero components of the TT tensor h 
are h^^,, where fi,v — t,x,y^z. This gauge is also called the Randall- Sundrum gauge. The advantage of this gauge 
(in contrast to the TT conditions applied to the 'complete oscillations', il^hab), as we will see below, is that for the 
graviton mode h^^, is independent of u and v. Then, at least for a > the conformal factor fl^ insures localization. 
In general, gravitational excitations are dependent on these variables, as well. They may or may not be localized to 
the intersection of the branes. In what follows, we will drop the superscript from the background metric, "^ab- 
The components of h satisfy the following wave equation 



V^VAh+-Ah = 0, 



(32) 



where A = t,x, ...,u,v and V"^ is a covariant derivative in the background metric, jab- Equation p2p is separable. 
Multiplying by the conformal factor we obtain 



9- 



1. 



2 V,,V,/i+-V""V„./i+-A/i = 0, 



(33) 



where the covariant derivatives are taken in the four dimensional metric g^i, and Ui is a shorthand notation for u 
and V. 

Now the first term of ()33|) contains only derivatives and metric components dependent on t, x, y, and z, while 
the rest of the terms depend on u and v only. Using a factorized ansatz for the wave function and introducing the 
gravitational mass m we obtain the following two equations (Here we assume a specific form for the 4 dimensional 



metric: Minkowski with a warp factor gyy 
a spatial coordinate, for dS4) 



Ox 



-gtt 



-2bz 



for AdS4 and gi; 



-2bt 



, where i denotes 



m h 



rn?h 



Aafl 



in 



dyg{u,v) 



+ 



4a 
dlh 



8a Ag(u,v) 
+ {u ^ v). 



4g(w,M) 9{u,v) \g{u,v) g{v,u) 



1 



1 



duh 



(34) 



(35) 



where the upper sign corresponds to dS4 and the lower one to AdS4. We use ^ = z or t and A = ±36^. 
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The solution of ((34l) is straightforward. It describes wave oscillations 



where kt = y/c^ + fc^ + fc^ + m^, provided the coordinates are much smaller than the curvature radius, 1/6, otherwise 

the wave is modified, as we will discuss it below. In particular, the m = mode always exists since /i2=const. is a 
solution of (|35[) . It describes a graviton, as long as the corresponding wave function is normalizable. In what follows, 
we will examine the solutions of (I35|) and their normalization. 

The exact solution of jMl) for dS4, using an ansatz h = e'(^'=-+y'=B+2fe=)y[t]^ is 



/[i] =e^^*/*i7i2)(e''*fc/5) (37) 



where Hn'^ = J„ + iYn is the Hankel function of the second kind, n — i^Jm? /IP- + 31/16 and k = \/k^ + k^ + k'^ 



When one uses the asymptotic expansion for large order but fixed argument to order (i.e. taking the limit b —>■ 0) [15| 
then we can write, after factoring out a diverging constant 

f[t] ^ e-''^^^+^+'^^''/''l (38) 



As a final note, as we indicated earlier, we mention that (|35p illuminates the special role of the seemingly arbitrarily 
chosen conformal factor, fl. The ground state wave function (m = solution of ([35]) ) is constant. Note however that 
the metric, and with that h was defined in ([T]) with a factor il^ extracted. Thus the true ground state wave function 
(provided it is normalizable) is exactly for every choice of our parameters. 

A. Mass spectrum of excitations 

In general, ([55]) is too complicated to be transformed to the volcano potential form, familiar in 5 dimensional brane 
theories, and solved using standard methods for solving a one dimensional Schrodinger equation. Consider, however 
that the coordinates a — u + v and S — u — v are somewhat analogous to the radial and angular coordinates in a 
Schrodinger equation. Thus, one expects, that solutions dependent on the 'radial coordinate', a, only well represent 
the spectrum of excitations. Among others, the gravition state, with a constant wave function, is clearly such a state. 

Indeed, using ([3]) we can see that for functions, independent of S ([35|) reduces to 

- Sh" + sign(a) ( ^X^-^ - -^h^ - aQs] h' = m^h, (39) 
\2a 4a / 

where function S is given in Table 1. Notice that the coefficients of pQ]) are independent of S, just like in a Schrodinger 
equation with rotational invariant potential the equation for rotational invariant states is independent of the angles. 
The reason for this cancellation is that the explicit i5-dependence of ([55]) comes only through multipliers A{6, a), which 
in turn are associated with a differentiation with respect to S. 

Before we can bring ((39)) to a Schrodinger form we need to introduce a new variable, uj, defined by 



UJ = sign(cr) / 
Jo 



l"l ds 



The integration in ([40|l can be analitically performed in the case of Solution I only. After the change of variables we 
can eliminate the first derivative from (j39p and obtain a Schrodinger equation with a volcano potential. This is done 
by factoring out an appropriate function p{uj) from h. Finally, we obtain the eigenvalue equation with the volcano 
potential in the form 

- d^h + V(uj)h - S(u;)^\2\ - Ag - + 2aS'mh = m^h, (41) 
4a 

where 

We use the notation h for the eigenfunction, though strictly speaking the function h in (|¥T|l differs by a factor of p 
from that used in 
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In and in S the variable a is understood to be substituted by a{u}), the solution of ([40]) . The explicit form to 
the potential in terms of variable lo can only be given for Solution I, but using (|4ip one can qualitatively analyze the 
spectrum for Solution II, as well. 

The range of variable uj and the behavior of the potential at the end of the range, loq, plays an essential role in the 
nature of the mass spectrum. Naturally, that range is dependent on the choice of gauge, but the variable lo is the 
natural coordinate. Anyway, the mass spectrum itself is gauge invariant. 

The m = eigenfunction of (|39p . candidate for the graviton as long as it is normalizable, is a constant. Then the 
corresponding eigenfunction of (j4ip is proportional to and given by 



ho{Lj) = 




where J7 = 1/w, S, and S' must be regarded as functions of cu. Here S' is the derivative of S with respect to aa. 
Note now that ho has no nodes, so provided it is normalizable it is the graviton state and the ground state of the 
eigenvalue equation. 

One can make a general comment concerning the mass spectrum. The potential depends on A only through the 
dimensionless parameter r that characterizes the function iS'(cr) completely. The remaining parameter a just sets the 
scale for the excitation spcctum. The value of r fixes the asymptotic curvature of the six and four dimensional spaces. 
Introducing the parameter ro = — 1 for Solution I and rg ~ —50/33 for Solution II we have 

• dSg and dS4 for r < tq 

• Mg and dS4 for r — tq 

• AdSe and dS4 for > r > tq 

• AdSe and M4 for r = 

• AdSe and AdS4 for r > 

B. Normalization 

The normalization condition is very simple in the Schrodinger type, 'w-representation' .The norm is finite if the 
integral over the square of the wave function converges. It is sufficient to know the wave function at the endpoint of 
the range of w, too (possibly infinity), as only the singular behavior at ljq can make the solution non-normalizable. 
Even when we are not able to solve the eigenvalue equation exactly we are able to determine the behavior of the 
solutions at ujq. However, the true restriction on the eigenfunctions is stronger. We must consider only solutions 
that can be orthonormalized, otherwise the small oscillation cannot be expanded into modes of definite mass. The 
condition of orthogonalization is that the Wronskian, ^11/12' — /i2^i' vanishes at w = wq If the combination hih'2 
is infinite then this condition is certainly violated. This condition is particularly important in the case of inverse 
square potentials, which will often make their appearance below. 

Every case in which the range of ui is finite the potential energy has the endpoint behavior V ~ —j{ujo — (^)~'^- 
The only acceptable solution behaves (neglecting subleading terms of the potential) as a Bessel function, ^/xJQ{mx), 
near x = ljq — lj = Q. The Wronskian vanishes for a pair of such solutions at different values of m and the solutions 
are orthogonal. The Wronskian for two singular solutions, which behave like ^/xY{){mix) , diverges logarithmically at 
x = making these solutions unacceptable. Having a unique acceptable solution implies that imposing the boundary 
condition aX, lo — Q introduces the usual quantization condition and the spectrum is discrete. Then, owing to the fact 
that in terms of the original variables, u and v a zero mass solution (independent of u and v) exists, there is a physical 
graviton state. 

Note that the normalizability and orthogonalizability of the solutions of the Schrodinger-like equation in variable lo 
does not insure the finiteness of the contribution of the corresponding small oscillations to quantities like the 'energy' 
defined by Mannheim (l6j . That requires the convergence of integrals over the variables u and v of the covariant 
weight, and the square of the small oscillation. This constraint come into focus for a < solutions, for which, 
depending on r, the end of range of w = 1 + a(\u\ + \v\) may he w = = 0, at which point j dS^fg ~ has a 
singularity. 

In tables 2. and 3. we summarized our results for both of our solutions for a > 0. The nature of the solutions 
changes as a function of r = — A/a^, among others the curvature of the 6 and 4 dimensional spaces, the question 
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TABLE II: Solution L Curvature of space, range of u), properies of the excitation spectrum and the range of weakly coupled 
states, as a function of parameter r. 



r 




space 


Graviton state Level spacing 


Weakly coupled 


r 


> 


AdSe, AdS4 


isolated 


Am ~ (-A)i/2 


m < (-A)i/2 


r 


= 


AdSg, M4 


marginal 


continuous spectrum 


m < (-A)i/2 




1 <r 


< AdSe, dS4 


isolated 


Am - (-A)i''2 


m < (-A)i''2 


r 


= -1 


Me, dS4 


isolated 


Am - (-A)i''2 


none 


< 


= -1 


dSe, dS4 


isolated 


Am - (-A)i/2 


none 



TABLE III: Solution II: Curvature of space, range of uj, properties of the excitation spectrum and the range of weakly coupled 
states, as a function of parameter r. 



r 


space 


Graviton state Level spacing 


Weakly coupled 


r > ro 


AdSe, AdS4 


isolated 


Am - (-A)i/2 


none 


r = ro 


AdSe, AdS4 


isolated 


continuous spectrum none 


ro > r > 


AdSe, AdS4 


isolated 


Am - (-A)'/2 


m < (-A)i/2 


r = 


AdSe, M4 


marginal 


continuous spectrum 


m < (-A)i/2 


-50/33 < r 


< AdSe, dS4 


isolated 


Am - (-A)i/2 


m < (-A)i/2 


r = -50/33 


Me, dS4 


isolated 


Am - (-A)i/2 


none 


-50/33 


dSe, dS4 


isolated 


Am - (-A)^/2 


none 



whether the graviton state is isolated, the order of level spacing of the spectrum, and the range of mass excitations 
with a weak coupling is displayed. The last two columns are important only for A << A, i.e. |r| << 1. This is the 
range in which corrections to the short and long range behavior of gravity do not contradict to observations in our 
current universe. 



C. Mass spectrum for a < 

If S{a) has no zero in the interval < w < 1 (here w — 1 — \a\a) then the solution is always unacceptable due 
to the divergence of the measure at w = 0. In all cases when a zero exists the parameter r is large. Though these 
solutions are mathematically acceptable, they could not be used for phenomenological purposes, for reasons discussed 
in the a > case. 



IV. CONCLUSIONS 

We solved the six dimensional Einstein equations and Israel junction conditions in a space with mirror symmetry 
and exchange symmetry, (.^2)^, in the two extra dimensions, when the 4-dimensional space at the intersection of two 
5-dimensional branes is warped. Due to their coordinate dependence, the junction conditions can only be solved if 
the gauges in the six and five dimensional actions are misaligned. We found two solutions, one with induced gravity 
terms on the branes, and one with an induced cosmological constant only. For each of these solutions we found two 
possible alternatives, having conformal factors that are increasing and decreasing, respectively, as the modulus of the 
5th or 6th coordinate increases. 

We also studied small perturbations around our solutions. Due to the complexity of the equations satisfied by 
these perturbations we restricted ourselves to those depending only on the sum of the extra coordinates. One of the 
oscillatory solutions of these equations is the massless graviton state, which is always localized at the intersection of 
the 4-branes. In addition there is an infinite set of solutions corresponding to gravitational excitations. The analytic 
form of the 'volcano potential' could only be derived for the solution with induced 5 dimensional gravity, but the 
graviton state and the nature of the gravitational excitation spectrum could always be analyzed. 

We found that for increasing conformal factors the graviton state is not normalizable. This is similar to the situation 
found in 5 dimensions by Karch and Randall [6]. However, here the lowest lying normalizable state is massive, leading 
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to an unacceptable short range gravitational force. 

For decreasing conformal factor, depending on the values of the appropriately normalized 6-dimensional and 4- 
dimensional cosmological constants, Ag and A, we obtained a variety of spectra. For most values of the cosmological 
constants we found an acceptable spectrum, with a zero mass graviton localized at the intersection of the two branes 
and a gravitational excitation spectrum that can be tuned not to lead to contradictions with experiments verify- 
ing Newton's law in the laboratory, provided Aq >> A. This work is an extension of similar investigations in 
5-dimension Q Q to 6 dimension. Previous studies of a 6 dimensional model with similar symmetries allowed 
only for M4 space at the intersection of the branes 

The model we studied requires the input of the value of the 4- dimensional cosmological constant by hand. We find 
this esthetically not very satisfactory. This blemish could in principle be circumvented by the addition of a global 
scalar, or several scalars, like quintessence, which could possibly generate a cosmological constant which may even 
depend on time. We intend to return to this problem in a future publication. 

We are somewhat perplexed by the existence of two acceptable solutions at given values of the input parameters 
(cosmological constants). We have not been able to find a physical principle which would allow us to discard one of 
the solutions yet. One possibility is to investigate observable post-Newtonian effects. The models under discussion 
have two separate gravity actions in 6 and 5 dimensions. Therefore, jus t like the Dvali-Gabadadze-Porrati (DGP) 
model [12], they could suffer from the van Dam-Veltman-Zakharov [l3l. fl3| ty pe of discontinuity and could be severely 
constrained. For more details we refer the reader to the following papers [l7l[l8l[T9l.[20j. 

The research reported here was supported in part by the U.S. Department of Energy grant No DE-FG02-84ER40153. 



APPENDIX A: SOLUTION OF THE GLOBAL EQUATIONS 



The two Einstein equations that we use to find the solution are then 

' lOa^ lOa^ 



Gun = 9uu -7 + T - 2X{1 + a{u + v)) 

g{u,v) g{v,u) 

2a (g^u{v,u) g,v{v,u) 



{l + a{u + v)) 



9iv,u) V giu,v) g{v,u) 
2a ( g,viu,v) g,u{v,u) 



\ + aiu -\- v) \ g{u,v) g{v,u) 

g ,u{v.u 
g(v,u) 

1 1 



0, 



(Al) 
(A2) 



where A = 87rG4A. Substituting ^ffjl^ from into ([JTI 

we obtain a first order differential equation for the 

symmetric function 



F{u,v) 



'^diUyV) 25(1), U) 



(A3) 



20a^F{u,v) - 2X{1 + aiu + v)f + Ae - {I + a{u + v))4:aF^{u,v) = 



(A4) 



First of all, this equations implies that the derivative F^y{u, v) is also symmetric. That is possible only if F{u, v) is a 
function of ct = u + f only. Keeping this in mind it has an analytic solution 



F{u,v) = S{a) = 4 
a^ 



il + aaf^^^ + C{l + aaf 



(A5) 



where C is an integration constant. 

Now defining S = u — v and using (jA3[) we can write g{u, v) in terms of an odd function of S, A{S, a) as 



g{u,v) 



S{a)+A{S,a) 

Inserting (jA6| into (jA2| we obtain a differential equation for A{S,a) 

A{S, a)A,,iS, a) + S{a)[A^s{S, a) - S'{a)] = 0, 



(A6) 



(A7) 
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We have not been able to give a solution of (|A7p in closed form, but the terms of its Taylor series in the variable 5 
can be readily calculated. We will write 

oo 
fe=0 

Then by substituting (jASp into (|A7|) we obtain 

ao{a) = S'ia). (A9) 

and the recursion relation 

^ k-l 

""'^"^ = - (2k + l)S(a) ^ (AlO) 

Using the recursion relation we can easily generate the coefficients in arbitrary order. We obtain 

S'ia) S" (a) 



ai{a) 



3S{(j) ■ 

S'{a) (S'iaf S"{a) + 2 S{a) S"(af + S{a) S'{a) 



=- ,,fy,5 ^S'{afS"{a) + nS{afS"{af-%S{a)S'{afS^^\a) 
+ 2&S{af S'{a)S"{a)S^^\(T) + S{a)S'{<Tf (--2QS"{(jf + iS{a)S^^\a) 



For the special case of Ae = 0, i.e. assuming that the six dimensional cosmological constant vanishes, still we get 
a curved space on the 3+1 dimensional word. In fact, at least in solution 1, A[6^a) simplifies considerably. We will 
consider this case next. 

First of all, S{cf) = w^, where we introduced the notation w = 1 + aa. If we also introduce the notation 

(AH) 

and we regard A as a function of w and z, then A can be written in the form 

A^u?h{z), (A12) 

where h{z) satisfies the ordinary differential equation 

2h^ + hh'{l -z)-h' = 2 (A13) 

We know the expansion of h{z) around z — \ already. It is quite interesting to investigate h{z) around 2; = 0, as well. 
The exact point, z = is not physical as a(5 < 1 + aa, but at large u and small v (near the u axis) z approaches 0. 
Numerical solution of this equations shows that /i(0)=l and h'{l) = 0. Of course we know that /i(l)=0. Then we can 
write h{z) = 1 — k{z) and expand (|A13p in k, keeping linear terms only. Simple calculation shows that (|A13p reduces 
to — 4A: + zk' = leading to h{z) = 1 — cz'^ + O(z^) (here we also calculated the next to leading term of k). This form 
is consistent with the expansion of (|A13p in a power series of k. The constant c can be identified from the numerical 
calculation as c = 1/4. 

The case Ag 7^ can also be investigated analytically near z — Q (near the axes, far from the origin). Using the 
ansatz 

A{5, a) = S{w) - h{z)f{w) (A14) 
(jA7p can be expanded to retain the first power of h only. One obtains a separable equation 

h{z)[S{w)f'iw) + f{w)S'{w)] - S{w)f{w)h'{z)z = (A15) 

The solution of this equation is 

S{w) 

h{z) = cz". (A16) 
Note that in the special case S = uP' the integration constant n — 4. 
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APPENDIX B: JUNCTION CONDITIONS WITH GAUGE TRANSFORMATIONS 



When we perform (jl2p the components of the metric tensor, g^^^i Qw and g^v change in the foUowing way (the 
time-Hke and ordinary space-hke components are unchanged aside from the replacements of the variables ) 

5™ T"^ — T^- — 9uu ^ R[g{u,v)[h^y{v,u)f + g{v,u)[h^y,{v,u)f] , 
1 + a[u + vjY 

9vv = TT- — 7^, — 77:^ ^ 9vv ^ R[g{v,u)[h,u{u,v)]'^ + g{u,v)[h^„{u,v)]'^] , 
guv = ^ guv = -R[g{v,u)h^u{u,v)h^uiu,v) + g{u,v)h^^{v,u)h^y{u,v)], (Bl) 
where we use the notation g{u, v) — g{h{u, w), h{v, u)), g{v, u) — g{h{v, u), h{u, v)) and where 

}_ (B2) 

[h uiu, v)h^y{v, u) — h y{u, v)h u{v, m)]^{1 + a[h{u, v) + h{v, u)]}^ 

Now notice that the components of the gauge transformed metric tensor at u = w = are given by ^uu (0,0) = 
(jvviO, 0) = 5(0, 0)[h^u{u, v)\u=v=o]'^ ■ Since normalization requires that these quantities are equal to 1 we have imposed 
the condition h^u{u,v)\u=v=o = 1 on (fT2|) . 

In this Appendix we will prove that the gauge transformed junction contributions of the six dimensional theory 
can be cancelled by the five dimensional curvature scalar and cosmological constant terms. First we need to find the 
gauge transformed form of the junction contributions ([9]). We will concentrate only on the contribution at the v = 
brane as once we solved the junction conditions at the v = brane the junction conditions at the u = brane are 
automatically satisfied considering the symmetry of gauge transformations (fT^ . As a reminder, we must take linear 

combination of the ARuu, ARvJ and ARuv contributions, such that they correspond to functional derivatives with 
respect to the untransformed metric components. These combinations are 

A 7,(6) 1 + auf 3 8ag{h, 0) - (1 + a h)g^y{h,v)\y^o 

A„4l) = -^""^TT^^(-)[^-(-'")l-"]'^^f^ (B3) 

where we used the abbreviated notation h = h{u, 0). 

Now we are ready to substitute pip and (jB3p into junction condition (fT5|) . Note that the two conditions {Ji{u) — 
and J2(w) = 0) at AB — ^fi and at AB — uu contain two independent functions: k{u) = h{u,Q) and f{u) — 
h^v{v,u)\y=o. These two functions are unconstrained aside from the following conditions: fc(0) = 0, fc'(O) = 1, 
/(O) = 1. /(u) appears in a very simple manner in the two junction conditions as shown by (|B3p . We can easily 
eliminate this function from the two junction conditions by taking the ratio of the two junction conditions. The 
resulting condition (call it F{k{u),u) = 0) is a transcendental function of k{u) (not a differential equation) and, in 
principle, can be solved for it. It has the form 



'>ag{k,0) ~ {l + ak)g^y{k,v)\y=o _ + o-u) - A5 - 2[g[«,o)]2 



8agik,0) 2X{1 + au)^ ~ St: G 5 A5 



(B4) 



Once we solved for k{u) we can calculate f{u) as well. All this can be done provided the three constraints on k{u) 
and f{u) are satisfied. Note that due to these three constraints the junction conditions Ji(0) = 0, J2(0) = 0, and 
dF{h{u),u)/du\u=o = become independent of the form of the gauge transformations. All other (higher) derivatives 
of the junction conditions contain unconstrained derivatives of the gauge transformation function. These provide 
equations for those unconstrained derivatives. We can conclude that the junction conditions can be solved for the 
gauge transformations provided three constraints among the constants of the theory are satisfied. The conditions 
Ji(0) = and J2(0) — were given earlier in ([T4|) and (jl5p . while the condition dF{fi{u),u) /du\u=o — was given 
in im. 
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